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Abstract 

Many models of supersymmetry breaking, in the context of either super- 
gravity or superstring theories, predict the presence of particles with Planck- 
suppressed couplings and masses around the weak scale. These particles are 
generically called moduli. The excessive production of moduli in the early 
Universe jeopardizes the successful predictions of nucleosynthesis. In this 
paper we show that the efficient generation of these dangerous relics is an 
unescapable consequence of a wide variety of inflationary models which have 
a preheating stage. Moduli are generated as coherent states in a novel way 
which differs from the usual production mechanism during parametric res- 
onance. The corresponding limits on the reheating temperature are often 
very tight and more severe than the bound of 10^ GeV coming from the 
production of moduli via thermal scatterings during reheating. 
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1 Introduction 



In supergravity and (super)string models there exists a plethora of scalar and fermionic 
fields (we will loosely call them moduli $) with masses of the order of the weak scale and 
gravitational-strength couplings to ordinary matter. If produced in the early Universe, 
such quanta behave like nonrelativistic matter and decay at very late times, eventually 
dominating the energy of the Universe until it is too late for nucleosynthesis to occur. 
This poses a serious cosmological problem Typical examples of these dangerous 
relics are the spin-3/2 gravitino (the supersymmetric partner of the graviton), the scalar 
moduli which parametrize supersymmetric flat directions in moduli space and seem 
almost ubiquitous in string theory, and their fermionic superpartners. 

The slow decay rate of the <I>-particles is the essential source of cosmological prob- 
lems because the decay products of these relics will destroy the ^He and D nuclei by 
photodissociation, and thus successful nucleosynthesis predictions p, §]. The most strin- 
gent bound comes from the resulting overproduction of D + ^He; this requires that the 
^-abundance relative to the entropy density s at the time of reheating after inflation 
§ should satisfy § 

— < io-'\ (1) 

s 

where the exact bound depends upon the mass m$. 

The dangerous moduli can be produced in the early Universe in a number of waysQ. 
They may be generated by thermal scatterings in the plasma during the process of 
reheating after inflation when the energy density of the Universe gets converted into 
radiation. The number density n$ at thermalization is estimated to be of the order 
of ~ 10~^ (Tr/Mpi), where is the flnal reheating temperature and Mpi = 

1.2 X 10^^ GeV is the Planck mass. This gives a stringent upper bound on the reheating 
temperature after inflation of about 10^ GeV 

The generation of the dangerous relics may also proceed through non-thermal effects. 

Scalar moduli fields may be copiously created by the classical gravitational effects on the 

vacuum state between the end of inflation and the beginning of the matter/radiation 

phase 0, ^ ||. In such a case, the upper bound on the reheating temperature is 

very dependent on the inflationary model, but it may be as low as 100 GeV. Another 

possibility is represented by the parametric excitation of dangerous relics occuring after 

inflation because of the rapid oscillations of the inflaton fleld(s). Gravitino production 

is an interesting example of this kind ^, 10, 11] since the non-thermal generation of 

^Hcre we will study only production of moduli quanta and we do not consider the possibility of 
moduli oscillations around their minimum. 
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gravitinos in the early Universe can be extremely efficient and overcome the thermal 
production by several orders of magnitude^. 

The scope of this paper is to show that the abundant production of dangerous relics 
is an unavoidable consequence of a period of preheating |]13[ after inflation. The corre- 
sponding limits on the reheating temperature are often very tight and more severe than 
the bound of 10^ GeV coming from the production of moduli via thermal scatterings 
during reheating. 

During the ffist stage of preheating nonperturbative quantum effects lead to an 
extremely effective dissipational dynamics and explosive particle production. Particles 
coupled to the infiaton field can be created in a broad parametric resonance with a 
fraction of the energy stored in the form of coherent infiaton oscillations at the end of 
inflation released after a few oscillations. This rapid transfer of energy density populates 
the Universe with quanta in highly nonthermal states which can be viewed either as 



classical waves travelling through the Universe |14| or as quantum particles in states 



with large occupation numbers. In the second stage of preheating, called semiclassical 



thermalization |14|, rescatterings of the produced fluctuations smear out the resonance 



peaks in the power spectra and lead to a slowly evolving state in which the power spectra 



are smooth |T^, |T5|, |T6[. The system begins to exibit a chaotic behaviour characteristic 
of a classical non-linear system with many degrees of freedom. In the course of the 
subsequent slow evolution, the power spectra propagate to larger momenta, eventually 
leading to a fully thermalized state. 

In supersymmetric and superstring theories coupling constants and masses appearing 
in the Lagrangian have to be thought as functions of the dimensionless ratio $/Mpi. 
For instance, a generic coupling constant h is in fact a function of the scalar moduli 

(2) 

where c is a coefficient usually of order unity and 5$ = $ — ($). This expansion 
introduces a gravitational coupling between the scalar moduli and the matter fields. 
If the underlying theory is supersymmetric, similar couplings of the scalar field $ and 




^ Recently, this statement has been criticized in ref. ||T^ , where it was shown that, in a simple Polonyi 
model, gravitino production is suppressed. Undoubtedly, as emphasized in refs. 0, ^, |l0[ 0, the non- 
thermal gravitino production depends sensitively on the assumptions on the inflationary model and its 
couplings to the supersymmetry-breaking sector. Non-thermal production of gravitinos is expected to 
be significant and to pose a cosmological problem if - after inflation - there is a strong mixing between 
the infiaton sector and the supersymmetry-breaking sector (of the present vacuum). As it will become 
clear in the present paper, this is often the case because, during the first few infiaton oscillations when 
the value of the infiaton field is close to Mpi, generic Planck mass-suppressed interactions can give a 
large contribution to such mixing. 
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its fermionic superpartner to the matter fields will arise in the superpotential or in the 
Kahler potential. 

Since moduli fields couple to any form of matter at least gravitationally, the produc- 
tion is present in a wide variety of inflationary models which have a preheating stage. 
In general, the equation of motion for the (quantum) fleld 5$ during preheating is 

□5$ + 5$ = ^ h,Ji , (3) 

i 

where the currents Jj are functions of the (classical) fields created at preheating and 
have interaction vertices parametrized by the coupling constants hi. In a cosmological 
setting particle creation in a time-varying classical background is usually caused by the 
time dependence of the effective mass of the field (consider, for instance, parametric 
resonance or gravitational creation of particles). As a result, particles are created in a 
squeezed state. On the contrary, creation of moduli is caused by the non-zero classical 
current and therefore they are created as coherent states. From a quantum mechanical 
point of view, this production mechanism is novel and differs from the more traditionally 
studied generation mechanism of squeezed states. Furthermore, the moduli production 
does not depend upon the particular properties of the modulus unlike the gravitino 
production [|^ ^ |10|, [l^] and it is true both for scalar and fermionic moduli. 



In this paper we perform a calculation of moduli production during the preheating 
stage in various inflationary models. In sect. 2 we consider the cases of chaotic inflation 
with quartic and quadratic inflaton potentials, with and without couplings of the inflaton 
to a new scalar fleld X. The case of hybrid inflation is discussed in sect. 3, while sect. 4 
contains a study of the production at preheating of modulinos, light fermions with only 
gravitational-strength couplings. In sect. 5 we summarize our results. 



2 Moduli Production at Preheating in Chaotic In- 
flation 

Let us flrst consider a simple model of chaotic inflation with quartic potential V{(j)) = 
^(j)"^. As discussed in sect. 1, we expect a generic coupling of the scalar moduli quanta 
5$ with the inflaton fleld of the form AJ5$, where J = ccf)^/ (4Mpi). Let us denote the 
amplitude of the oscillating inflaton fleld as 00, while the frequency of the oscillations 
is ~ \/A0O; where 0o is the zero-momentum mode of 0. 

In this model fluctuations of the inflaton field do not grow during chaotization (they 
actually decrease due to the redshift) so the most efficient production of moduli fields 
take place at the end of the resonance phase. As we noted, one can describe the states 
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we are considering either as a collection of interacting classical waves or as a collection of 
particles in modes with large occupation numbers. Therefore at late times it is possible 
to treat the system classically and to solve the relevant equations of motion on a lattice. 
The initial conditions for the classical problem are found solving the preceding quantum 
evolution by means of the appropriate Bogolyubov transformations (for more details see 
Ref. 0). 



The results of the numerical integration on a 128'^ lattice for A = 10^^^ are summa- 
rized in Fig. 1, where we plot the values of (0^), ($^) and of the inflaton zero mode 
squared 0q as a function of the conformal time and in units M|j. The field variance 
is given by (0^) = [(27r)'^V]~"'^ / (i^k0^, where V is the volume and 0k is the Fourier 
k-mode of the field 0. Here and in the following we have taken c = 1. Time and parti- 
cle momenta are given in units of \/X(f)o{0), where 0o(O) is the initial value of the zero 
inflaton mode at the end of inflation. The time t = 0, when inflation ends, is chosen 
as the moment of the first extremum of 0(t). The maximum value of (0^) is achieved 
at the time tr, which marks the end of the resonance phase, and, as shown in Fig. 1 is 
about 10-^ M|i, corresponding to ($2) ~ 10"^^ M^j. At the same time the zero mode 
is about 4>o{tr) ~ 10"^ Mpi. 

We will now present some analytical estimates which approximately reproduce the 
numerical results of Fig. 1. To obtain the intensity of the moduli gravitational pro- 
duction which accompanies the creation and annihilation of the fluctuations generated 
during preheating, we first express the solution of the modulus equation of motion in 
terms of the retarded Green function 

= $i„(a;) + A / d'^x'Gretix - x')J{x'). (4) 



Here is the initial configuration at time t = and the retarded Green function is 
given by 

G„,(x - X') = 0(t - t') I -^ """"''-''' e-^-' (5) 

In the case under consideration, the source term is given by J = c0^/(4Mpi). Expanding 
the fields in Fourier modes, performing the integral over spatial coordinates, and taking 
00 ^ 50k (justified by the numerical results of Fig. 1), we obtain 

Mpi Jo Uk 

A crucial point is that the effect is linear in the inflaton fluctuations 50k. Therefore, 
the phenomenon of moduli production is due to the oscillations (similar to oscillations 
of neutrino species in a medium) of the inflaton quanta into moduli in presence of 



4 



10-^ 
10-6 
10-^ 
10-8 
10-9 



10- 


10 




10 


11 




10- 


12 




10- 


13 




10- 


14 




10- 


15 




10- 


16 




10- 


17 





' I I I ij ' 



I [ 1 [ 1 [ 1 t I I [ 1 I I [ 1 T [ I [ 1 I I 1 [ 1 [ 1 t 1 n 




200 400 



600 800 1000 1200 1400 
T 



Figure 1: The variances of the inflaton and the moduli fields in units of M|i in a model 
with infiaton potential V = ^(p^ and moduh interaction cA0'^5<l'/(4Mpi), with A = 10"^^ 
and c = 1. The dashed line shows the amphtude square of the infiaton field oscillations 
in the same units. 
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the inflaton background (po, rather than to real rescatterings. From a pure quantum 
mechanical point of view, scalar moduli are generated as coherent states and not as 
squeezed states as in the usual parametric resonance. 

The fluctuations 50k grow exponentially during the resonance phase as e^* where 
2^ ~ O.O7\/A0o(O) and 0o(O) ^ 0.3 Mpi is the initial value of the inflaton field at 
the end of inflation. Therefore the integral in eq. (|^) can be approximately evaluated 
by assuming that the exponential determines the full time dependence of the integrand. 
Then we obtain that the maximum value of the variance of the modulus field, achieved 
at the end of the resonance phase, is given by 

(^^) - "^Kf'^'iif - lOV^ (0^). (7) 



Here we have used oOkiU) ~ VX4>o{tr) as the typical frequency at the end of the resonance 
phase. Taking from Fig. 1 the values 0o ~ lO'^Mpi and (0^) ~ 10-^M|i at t = U, eq. (0) 
gives ($^) ~ c^lO~^^M|i, in very good agreement with the numerical result. Notice that 
the dependence on A dropped out from eq. (J^). Our numerical calculation confirms that 
the maximum of ($^) is independent of the coupling A and it scales quadratically with 
c. 

Since it is a good approximation to assume that the dominant particle production 
occurs around the time tr, we find 



~ LL!k{tr){^'^) and ~ uJk{tr){<l)'^)- 

Using eq. (0) we obtain 



10-3Mpi 



at the end of the resonance phase. 

The ratio of the number density of moduli in units of the entropy density s ~ p^/^ ~ 
A^/'^0o(tr) is given by 

For A = 10~^^ (as fixed by the COBE normalization), this results contradicts the bound 
(P by about six orders of magnitude (for c of order unity) when the energy density in 
the scalar field is transferred to the energy density of a hot gas of relativistic particles. 
Notice that the ratio ^ does not depend on the time of thermalization, because both 
n$ and V^^^ vary with the scale factor a as a~^. 

In this calculation (and throughout all this paper) we have made the implicit hy- 
pothesis that the moduli fields are light during the preheating stage. However, the 
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Figure 2: The variances of the inflaton, X and moduh fields in units of in a 
model with potential V = ^(f)^ + ^0^X^ and moduli interaction cA0'^5$/(4Mpi), with 
A = 10^^^, q = 5'^/(4A) = 25 and c = 1. The dashed fine shows the amplitude square of 
the inflaton field oscillations in the same units. 

same generic expansion in eq. (2) is expected to provide the moduli fields with a mass 
of the order of the Hubble rate H during the preheating stage. Consider the simple 
model of chaotic infiation with potential V{(j)) = ^(j)^. The expansion of the coupling 
A($/Mpi) at the quadratic order in $/Mpi leads to a term in the Lagrangian of the 
form ~ \{6^/Mpif(j)^ oc H^{6^f. Nevertheless, our numerical and analytical esti- 
mates are not affected by treating the moduli as light states. This is because moduli 
fields are generically produced with frequencies u much larger than the Hubble rate 
{lu ~ V^0o(^r) ^ H(tr) in the chaotic model A0^) and kinematically they may be re- 
garded as massless. This does not come as a surprise since it is known that superheavy 



particles may be easily generated at preheating |T7 . 

A simple variation of the model considered so far is the one in which a massless 
inflaton field is coupled to another field X, V = ^X4>^ + \g'^(l)^X'^. The oscillating field 
(f) amplifies fiuctuations of the field X via parametric resonance. The strength of the 
resonance can be parametrized by the g-parameter q = g'^/4:X. Fluctuations of the 
inflaton fleld grow rapidly because of the back-reaction and the energy of the inflaton 
zero mode is quickly released into quanta of the X-fleld and the inflaton fleld. 
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Figure 3: The variances of the inflaton, X and moduh fields in units of Mpj in a 
model with potential V = + ^^^X^ and moduli interaction g'^(j)^X'^6(^ / {2M-pi), 
with A = 10^^^, q = g"^ / (4A) = 25. The dashed line shows the amplitude square of the 
inflaton field oscillations in the same units. 

For q = 25, we numerically find that the maximum value of (0^) is about 5 x 10~^M|i 
while (X^) is about one order of magnitude smaller, see Fig. 0. The resonance stage 
ends when (plitr) ~ 5 x 10~^ Mp^ ^ (0^)- If "we again restrict ourselves to the interaction 
of the form cA5$0^/ (4Mpi) , moduli are produced through the scatterings of the infiaton 
quanta and - at the end of the resonance stage - ($^) ~ 10~^^ ^pi- 

The analytic estimate in eq. (|^) well reproduces this numerical result for ($^). The 
corresponding ratio derived from eq. (p!oD , violates the bound (1) by seven orders of 
magnitude for c = 1. 

Another option is to have a coupling of the form g'^6^(f)'^X'^/{2Mpi). In such a case 
the production of moduli cannot take place through oscillations at the initial stage of 
the parametric resonance, but only via inverse decays of already created fiuctuations of 
the X field. This is equivalent to say that the effect in the Green function method is 
not linear in the fluctuations, but at most quadratic. However, moduli are produced at 
the initial stage of the semiclassical thermalization and the result is not much different 
from the previous case, see Fig. |^. 

The main difference during the evolution arises because the effective couplings in 
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the system depend upon the avereage value of the fluctuations of the moduli field. This 
changes effectively the parameter q and the time development of the resonance which 
is very sensitive to any change of q. 

We can now use the Green-function method to give an estimate which approximately 
reproduces the results found in Fig. |. The Fourier modes of the moduli fields at time 
t can be written as 

^^^'^ = V&f '"^"'^'~'^^^ '^o(0 / d^vX,{t')X^-,{t'). (11) 

AIVlpi Jo UJ}^ J 

As expected, the moduli Fourier modes are quadratic in the fluctuations of the X 
fleld. We can perform the time integral by assuming that the leading time dependence 
comes from the exponential growth of Xp, which is proportional to e'^-^*, with — 
O.1a/A0o(O). Since the momentum distribution of the X fluctuations is sharply peaked, 
the integral over k of the square of eq. (^) at the end of the resonance gives 



9'<Plitr){X^) 



lovM^. (12) 



'PI 



Here we have used cuk ~ "\/A0o as the typical frequency. From Fig. ^we gather that, 
at the maximum, {X^) ^ lO-^M^i, corrisponding to ~ lO'^M^^. For q = 25, 

eq. ([I2|) gives ($^) ~ 10^^^M|i, in fair agreement with our numerical results. 

In all cases we also have calculated the particle occupation numbers to make sure that 
the lattice spacing and the integration box are suitably chosen and the relevant particle 
momenta are well represented. Particle number densities in the comoving reference 
frame are shown in Fig. ^ for the same model which is represented in Fig. ^ and at the 
time T = 180. Let us also note that the production of moduli during the preheating 
stage should be analyzed by letting the moduli be coupled to the matter flelds in all 
the interaction terms. Therefore, our estimates provide a lower limit on the abundance 
of moduli states. 

Reassured by the capability of reproducing the numerical results by analytical means, 
we have also considered the quadratic inflaton model with quartic coupling to another 
fleld X, V = |m^0^ + l^f^^^X^. The g-parameter is g = g'^(f)l{0)/4:m'^. It is also 
useful to introduce the redshifted resonance parameter at the end of the resonance 
stage qr = q{4>Q{tr) / 4>o{0)y ■ Parametric resonance can fully develop in an expanding 
Universe if ~ 1 UHl- Resonant production is most effective for fluctuations of flelds 
which couple not too weakly but also not too strongly, those with g,. ~ 1. For g^ 1, 



the maximal size of X fluctuations is suppressed by non-linear effects [0, |T6[. The zero 
mode rapidly drops at the end of the chaotization stage and all its energy is transferred 
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Figure 4: Particle number densities of X, inflaton and moduli fields, from top to bottom, 
at T = 180 in the comoving reference frame, in units of V^0o(O), for the same model 
studied in Fig. ^ 



to fluctuations the variance of the inflaton fleld is much larger than the value of 
the zero mode at the chaotization phase. 

We now want to perform an analytical estimate of the number density of moduli 
fields for moderate g^- If the mass term has a dependence on the modulus field 

= m^($/Mpi), this will introduce a trilinear coupling m^0^5$/(2Mpi) between the 
quanta of the modulus and the infiaton quanta. In such a case, the Green-function 
method previously illustrated gives a modulus variance at the end of resonance equal 

($^) Ml! (13) 



Here we have taken ~ m as the typical energy and /i ~ 0.2 m [T^. For ~ 1, which 
corresponds to g ~ 10^ ||15|, the infiaton completely decays into fiuctuations during the 



chaotization phase; at the end of it 0q ~ (0^) ~ n^/m ~ 10~^Mp[. Therefore, from 
eq. ([T3|), we obtain n^/n^ ~ 10^^. 

Notice that we could obtain the same estimate using the Boltzmann equation for 
the number density n$: 



171^ nl 



n^ + 3Hn^^—-^. (14) 
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Assuming again c<j ~ m and knowing that grows exponentially with rate 2/x, we get 

At the end of the reheating stage, the energy density stored in the infiaton is con- 
verted into a relativistic thermal bath with temperature ~ (mn^)^/^. The number 
density of moduli per entropy density will then be 

!^ ^ ^ io-^Il_ (15) 

s m m 

Since the COBE normalization gives m ~ 10^'^ GeV, the upper bound on the reheating 
temperature becomes 

Tr < 10^ GeV. (16) 

If the coupling of the modulus is 5'^5$0^X^/(2Mpi), the Green function method 
gives 

'9'<Pl{tr){X'y 



- 2 






0.2Mpi 



(17) 



Here we have assumed uo m and ~ 0.2m. Repeating the same argument used 
above, we find ~ 10''g^(X^)^Tr/(Mpim). Since at the end of the resonance phase 
Qr ~ 1, we find that the bound (1) gives a limit on the reheating temperature of 10^ GeV 
if - for instance - we use q = 10® and the maximum value of (X^) ~ lO^^Mpj [|l6l . 

So far we have always considered moduli couplings arising from the scalar potential, 
but similar results are also obtained when the moduli couplings are derived from non- 
minimal kinetic terms. For instance, from an interaction of the form 5$0^/Mpi we 
estimate 

In the case of chaotic inflation with quadratic potential, eq. (|T8|) is equivalent to eq. (|13D, 
valid for moduli coupled to the infiaton mass term. 

3 Moduli Production at Preheating in Hybrid In- 
flation 

In this section we study the moduli production during the preheating phase following a 
period of hybrid inflation . This scenario involves two scalar fields, the infiaton field 
0, and the symmetry breaking field a, and different mass scales and couplings. During 
inflation, the infiaton field (f) rolls down along a flat potential while the Higgs field a 
is stuck at the origin providing the vacuum energy density driving inflation. However, 
when (f) gets smaller than a critical value 0c, both fields roll down very quickly towards 
their present minima, completing the inflationary phase. 
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Figure 5: Variances of the inflaton, Higgs and moduli fields in units of f ^ in the hybrid 
model with /\ = 10, A = 10"*^ and v = 10~^Mpi. The moduli are coupled according 
to the g'^(j)^X^(^/{2Mpi) interaction. 

We take the simplest hybrid inflation potential as suggested by Linde [|18| f\ 

Vi<^, a) = ^ {v' - + lml<^^ + i/^/)^^ . (19) 

This potential has a valley of minima at a = for 0^ > 0^ = p-f ^. Most of the inflation 
occurs while is slowly rolling down from its initial value to 0c- During this phase 
(7 = and the potential reduces to ^ = ^v^ + |m|0^. For < 0c, the concavity of the 
potential in the a direction changes so that the 0-axis is a ridge of unstable maxima. 

Inflation ends because of the growth of quantum fluctuations in a. For > 0c, these 
are quickly damped out, but beneath 0c they cause the fields to fall very quickly to 
global minima at = and = f ^. There, they will execute damped oscillations. 

The results of numerical integration for the case v = 10~^Mpi are shown in Fig. ^, 
where we have supposed that the modulus is coupled to the inflaton field and the Higgs 
field via the coupling (yf^$0^(T^/(2Mpi). In hybrid models the coupling constant A does 
not have to obey the same restrictions as in the case of chaotic inflation. However, for 
the numerical integration we have chosen a rather low value, A = 10~®. This is because 
^For other hybrid inflation models, including those motivated by supersymmetry, see ref. B|. 
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the particle number in the inflaton or Higgs excitations at the beginning of chaotization 
is n{k) ~ 1/A and to be able to treat moduli fields classically, we need occupation 
numbers n$ ^ 1. Since n$ ^ nlk), we need a rather small value of A. For the run 
presented in Fig. |^ we have chosen g'^/X = 10. 

The initial conditions were chosen in such a way that Higgs oscillations occur only 
around one of the minima, and domain creation and other effects leading to complete 
chaotisation during the first oscillation (see Ref. can be disregarded. 

Our numerical findings show that the variance of the moduli fields does depend only 
upon the ratio g'^/X. 

The Green function method tells us that 

since the typical frequencies are given by uj"^ ~ X{a)^ and uj^. ~ g'^{a)'^. Using the same 
values of the couplings and mass scales adopted in Fig. 5 and reading from there the 
variance of the infiaton field (0^) ~ 10~^f^, we get ($^) ~ 10~^f^, in fair agreement 
with our numerical results. Notice that the numerical values of the infiaton and Higgs 
variances at their maximum are well reproduced by (0^) ~ {X/g'^)v'^ and (a^) ~ v"^. 
Using these values we estimate ($^) ~ {X/ g^){v^/Mpi), in good agreement with our 
numerical results and with the fact that these results depend only upon the ratio g'^/ X. 

If the final reheating takes place via decay of the lightest particles in the system, i.e. 
the Higgs particles with mass squared Xv"^, the final moduli number density-to-entropy- 
ratio is given by 



s ria y/Xi 



(21) 




T.<v/A r^ ^^^^M 10-GeV. (22) 



which implies 



This result illustrates how, in certain ranges of parameters, the coherent production of 
moduli can become more dangerous than the generation of moduli via thermal scatter- 
ings. 



4 Modulino Production at Preheating 

So far we have been considering the generation of scalar moduli during the preheating 
stage. In this section we consider the production of the fermionic partners ipq, of the 
moduli fields, sometimes dubbed modulini. As for the scalar moduli, the modulini are 
generically coupled to the matter fields through couplings suppressed by M-p\. 
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Consider - for instance - the superpotential W = -^0^ giving rise to the inflaton 
potential V = A0^. If the couphng A is a function of ^/Mpi, where now $ has to be 
considered as a superfield, the Lagrangian will contain the following couphng 

(23) 

Mpi 

where have denoted with the inflatino, the fermionic partner of the inflaton field. 

We can again make use of the Green function method to obtain the Fourier compo- 
nent of ip<s,. Using the retarded Green function for a fermionic field, we find 

ifj^Ck^t) = —— / dt' < -i7°sin[u;fc(t-t')] cos[uJk{t - t')]\ (f)l{t') ilj^{k,t'). 

Mpi JO I UJk I 

(24) 

We can then approximatly estimate the quantity {ip^4'<s>) as 

(^^^*) ^ ^4(^^^*)- (25) 

One can now use as a guide the recent results obtained in the theory of generation of 
Dirac fermions during and after inflation pO| . During the first inflaton oscillations, the 



Fermi distribution function of the inflatini is rapidly saturated up to some maximum 
value of the momentum k, i.e Uk — 1 for k < /Cmax and it is zero otherwise. The resulting 
number density is therefore ~ ^max- The value of /cmax is expected to be roughly 
of the order of the inverse of the time-scale needed for the change of the mass of the 
inflatini. In the model under consideration, such a mass changes by an amount V^0o 
in a time scale (-\/A0o)~^ and one expects /Cmax ~ VX(f)o- Therefore 

{^^^^) ~ X'/'<Pl{0), (26) 

where we have used the fact that most of the fermion production takes place during 
the first oscillation. Under these circumstances, one can also assume that the number 
density of inflatini n^^ and modulini n^^ is well measured by {ip^ip^) and {ip<s,'ip<s>), 
respectively. This gives 

0o(O) 

n^, ^ Ml, ' ^'^^ 
where we have taken u ~ fcmax in Eq- (11) • 

Finally, the ratio of the number density of modulini in units of the entropy density 
s ~ p^/^ ~ A^/^0g is given by 



W-^\ (28) 
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This result will contradict the bound (|1|) by about an order of magnitude when the 
energy density in the scalar field is transferred to the energy density of a hot gas of 
relativistic particles. 

Similar considerations apply to the infiaton model with quadratic potential. Sup- 
pose there is a coupling of the form m0^$/Mpi in the superpotential. This gives rise 
in the Lagrangian to a coupling of the modulini with the fermionic partners of the 
infiaton of the form Using the Green function method we again obtain 

^V* ~ '^V0('^o(O)/Mpi)2. This gives n^^/n^^ ~ (m/Mpi)^ ~ 10^^^ and the generation 
of dangerous relics is not large enough to provide a strong constraint on the reheating 
temperature. 

In the case of the hybrid model though, the bound is more restrictive. In super- 
symmetric hybrid models, the coupling constants A and g are related, and typically one 
finds A/(7^ ~ 1 0]. The COBE normalization implies that typically i; ^ 5 x 10^^ GeV. 
During the oscillations of the Higgs and the infiaton field, one expects that the number 
density of "Higgsini " ipcr and infiatini ip^ is of the order of the number density of Higgs 
particles excited during the preheating phase, n^^ ~ n^^ ^ ^ X^^'^v^. If the modulus 
couples to the Higgs superfield via the interaction \0\{(^/Mpi)a'^(f) in the superpotential, 
the number density of modulini is expected to be 



The corresponding bound on the reheating temperature can be obtained using the same 
arguments which led us to eq. (|2|); we find < (5 x lO^^GeV/v) 10^°GeV. 

5 Conclusions 

In this paper we have investigated the production of scalar and fermionic moduli fields 
which occurs during the preheating stage after infiation. If produced with too large 
abundances, the late decays of these relics may jeopardize the successful predictions of 
standard big-bang nucleosynthesis. Since the moduli couple to any form of matter at 
least gravitationally, their generation is unescapable in the large class of infiationary 
models which have a preheating stage. They are generated as coherent states initially. 
From this point of view, the production mechanism of moduli differs from the gener- 
ation mechanism of squeezed states during the parametric resonance which has been 
studied in recent years. We have shown that moduli are efficiently produced during 
the preheating stage. The corresponding upper bound on the reheating temperature is 
model dependent and sensitive to the details of the stage of thermalization following 




(29) 
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the period of parametric resonance. It is often tighter than the bound of 10^ GeV found 
considering the generation of moduh via thermal scatterings during reheating. 
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